Giant mesoscopic spin Hall effect on surface of topological insulator 
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We study mesoscopic spin Hall effect on the surface of topological insulator with a step-function 
potential. The giant spin polarization induced by a transverse electric current is derived analytically 
by using McMillan method in the ballistic transport limit, which oscillates across the potential 
boundary with no confinement from the potential barrier due to the Klein paradox, and should be 
observable in spin resolved scanning tunneling microscope. 



Topological insulator (TI) with time reversal invari- 
ance has recently been proposed theoretically and ob- 
served in experiments 1 - - — . In three spatial dimensional 
(3D) TI, electronic structure is characterized by a bulk 
gap and a gapless surface mode described by odd number 
of branches of Dirac particles, which is protected by time 
reversal symmetry. The surface states are helical, where 
spins are locked with momentum. The gapless Dirac dis- 
persion mode of the surface states has been confirmed in 
angle resolved photoemission spectroscope, and the pre- 
dicted spin-momentum correlation has also been reported 
in experiments 6 -^. Because of the strong correlation be- 
tween the spin and momentum, the surface states of 3D 
TI may be a potentially idea system to study spintron- 
ics, where the electron's spin degree is used to manipulate 
and to control mesoscopic electronic devices. It is inter- 
esting to note that the recent progress in TI is closely 
related to the development of the spin Hall effect (SHE) 
in the past several years, a sub-topic in spintronics. 

The SHE refers to a boundary (surface or edge) spin 
polarization when an electric current is flowing through 
the system. There have been extensive studies on 
the SHE in the conventional semiconductors or metals 
with spin-orbit coupling, both in experiment*^— and in 
theory 1 - 4 . - . 19 . The SHE is often classified into "extrinsic" 
(impurity driven) or "intrinsic" (band structure driven). 
Since arbitrarily weak disorder destroys the intrinsic SHE 
in 2D infinite system with linear spin-orbit couplin g 20 ' 21 , 
there have been considerable interests on the mesoscopic 
systems in the ballistic limit, where the disorder may 
be ignored 2 ^—. In the ballistic limit, the electric field 
is absent inside the system, and the spin polarization 
is resulted from spin precession around the lateral con- 
fined potential. The SHE of 2D semiconductor system in 
the ballistic limit has been studied theoreticall y 23 ' 24 ' 26 . 
While the ballistic spin accumulation is predicted near 
the potential barrier, the effect has not been observed 
in experiments for the weakness of the effect in realistic 
semiconductors or for the difficulties to detect the spatial 
distribution of the spin polarization in the sandwiched 
interface. The surface states of the TI represent a differ- 
ent type of 2D system where the spin-orbit coupling is 
strong, and the surface state can be probed directly by 
scanning tunneling microscope (STM). This may provide 
a new route in study of the SH E 28 ' 29 and spintronics in 
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FIG. 1: (Color online) Schematic illustration of the proposed 
ballistic spin Hall effect on the surface (x — y plane) of a 
topological insulator. A back gate of voltage V2 is applied 
at the right side with x > 0. Electric current flows along 
y-direction, resulting in spin polarization oscillation along re- 
direction across the barrier. 



general. 

In this work, we report theoretical prediction of a giant 
SHE on a surface of 3D TI with a step-function potential 
in the ballistic limit as schematically illustrated in Fig. 
1. By using McMillan Green function method, we de- 
rive analytic expressions for the electric current-induced 
spin polarization, which is oscillating across the poten- 
tial boundary, and is not confined by the potential due 
to the Klein paradox. For a typical TI, the amplitude 
of the local spin polarization is estimated to be as large 
as 20% at the Fermi level near the boundary, which is 
much larger than the SHE in a typical 2D semiconductor 
system, and should be observable in spin resolved STM 
experiments. 

We consider surface states of TI, described by an ef- 
fective Hamiltonian in the x-y plane, 



H = v F (p x er) • z + V(x), 



(1) 



where p is the electron momentum, <x are the Pauli matri- 
ces, and vf is the Fermi velocity. The system is transla- 
tional invariant along the y-axis, and has a step-function 
potential at x = 0, which separates two regions along 
the x-axis: region 1 at x < and region 2 at x > 0, as 
illustrated in Fig. 1, 
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Vi=0 if x < 0; 
V 2 if x > 0. 



where V2 is a constant. A voltage of V y is applied across 



2 



the surface to induce an electric current along the y- 
direction. We consider the ballistic limit, where the elec- 
tric field inside the surface is zero. We will first construct 
the retarded Green's function by using the scattering 
wavefunction, a method introduced by McMillan to study 
superconducting state. From the obtained Green func- 
tions, we calculate the local spin density in the presence 
of the electric current to show the profound SHE. The ex- 
perimental consequences and comparisons with the SHE 
in conventional 2D semiconductors will be discussed. 

The scattering wave functions can be constructed 
based on the eigen functions of the Dirac particle in 
Hamiltonian (1) in the two separate spatial regions. The 
eigen wavefucntions in region a corresponding to the en- 
ergy e and y-component momentum k y are given by 
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where %vpk ax — \/( e — ^a) 2 — (fivpk y ) 2 . By adjust- 
ing the gate potential V a relative to the Fermi energy 
Ep, the Dirac fermion carriers in region a can be tuned 
into electron-like (n-type, Ep > V a ) or hole-like (p-type, 
Ep < V a ) . Therefore the system may be viewed as n — n 
or n — p types of junction. 

The right (R) and left (L) moving scattering wavefunc- 
tions can then be found by using the standard transfer 
matrix method. For a n — n junction, we have 



<t>T(x,y) 

<t>T{x,y) 



/ <Pi(x,y) + r R n <p 1 (x,y) ifa;<0 

t^+(x,y) if.x>0. 

tiyi(x,y) ifa;<0 , 

¥>2 ( x >y) + r 2 n ft( x ^y) if»>0. 



where f^) L and T^) L are the transmission and reflection 
coefficients respectively, which are related by 
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where T is the transfer matrix, 
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The scattering wavefunctions for a n — p junction have 
similar form with that for the n — n junction, except that 

cf>2^ ■* are replaced by <p 2 ^ in the region of x > for the 
group velocity of a hole is opposite to that of an electron 
and that all the superindices of nn are replaced by np. 



We note that if k x is complex, the evanescent wave ap- 
pears. In this case, considering the asymptotic behavior 
of the evanescent wave, the scattering wave function for 
both n — n and n — p junction will have the form given 
in Eqn. ©. 

We are interested in the transverse effect of the charge 
and spin density profiles as an electric voltage is applied 
along the y-axis. To this end we construct the retarded 
Green's function, which satisfies the equation, 



(e - H)G r {x, x'; e, k y ) = S{x - x')I 



(4) 



where I is a 2 by 2 identity matrix. The solution for G r is 
a direct product^ of the scattering wavefunctions 4>r/l 
and the transposal wavefunctions (^l/ri 

c<<I)l{x, y)4> R {x',y) if x < x' 



c>0 fl (a;,y)^(x',2/) if x > x' 



(5) 



where c < and c > are the coefficients, which can be de- 
termined from Eqn. (4). Here, 4>l/r nas the same form 
as 4>l/Ri expect the replacement of the factor e lky ' y by 

g iky -y 

The local spin density of states for a giving k y (S) and 
the local charge density of states p at energy e can be 
found easily from G r , 



S(x;e,k y ) = ——lmTi[G r (x,x;e,k y )cr] 
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where the sum in p is over all the possible values of k y , 
and the contributions from the evanescent waves are also 
included. 

In the present case, the Green's function, hence the lo- 
cal charge and spin density of states can all be solved an- 
alytically. Here we shall focus on the spin z-component, 
which is most interesting and given by 



Tr[G r (x; e, k y )o z \ = F(e, k y )[e^* x Q(x) 
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where Q(x) is Heaviside function and S2 = sign(e — Vz). 
The charge density of states is directly related to the 
trace of G r , 
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FIG. 2: (Color online) Spatial distribution of local spin den- 
sity S z for a single Dirac fermion in Hamiltonian (1) with 
vf = 5 x 10 5 m/s. (a) n-n junction, Vi = 40 meV, and 
e = 60 meV. 9 = 63° (blue solid line) and 9 = 18° (red 
dashed line), (b) n-p junction, V2 = 60 meV and e = 35 meV. 
9 = 63° (blue solid line) and 9 = 36° (red dashed line). 9 is 
an incident angle, defined by e lB = -prw ■ 



The local spin density of states for a giving k v in the 
n — n and n — p junctions are plotted in Fig[2] One 
important feature is the non-confinement of the Dirac 
particle with the higher barrier potential (V2 > e) due 
to the Klein paradox. The local spin density strongly 
depends on the incident angle of the electron, as we can 
see from Fig. [2] The evanescent wave appears if |e— V%\ > 
hvFky > |e — V2I, and there is a critical incident angle for 
the condition of the evanescent wave. These features are 
typical characteristics of Dirac fermion, essentially the 
same as in the graphene. 

We now calculate the spin polarization along x- 
direction near the potential boundary x = induced by 
an electric current along the y-direction at zero temper- 
ature. We consider a voltage of V y /2 at the one edge 
and —Vy/2 at the other edge of the system along the y- 
axis, and consider the ballistic transport limit. The effect 
of the voltage at the two edges is to induce an imbal- 
ance of the occupied states between k y > and k y < 0. 
The states with k y > are occupied at energies below 
Ep + V y /2, and the states with k y < are occupied at 
energies below Ep — V y /2, with Ep the Fermi energy at 
V y = 0. By the time reversal symmetry, the local spin 
polarizations contributed from k y > and from k y < 
with the same energy cancel to each other. For small 
value of V y , we thus obtain the current-induced net spin 
density profile^^i^ 

ST(x; E F ) ~ \eV y \ £ S z (x; E F ,k y ) (10) 

ky>0 

To further analyze the current-induced spin polarization, 
we define local spin susceptibility Xz(x;Ep) and local 



spin polarization P z (x; Ep), 

Xz(x;Ep) = S?(x;E F )/^eV y 

P z {x-E F ) = Xz{x;Ef)/p(x;E f ). (11) 

The local spin polarization P z is a dimensionless parame- 
ter to measure the magnitude of the SHE. P z = 1 means 
the spins of the electrons at a space point x at the Fermi 
level are polarized. The experimentally measured local 
spin density is obtained by multiplying P z by the applied 
voltage and by the local density of states. In Fig|3] we 
plot Xz and the P z in both n — n and n—p junctions. The 
key features of the SHE in the system are summarized 
below. 1). There is a pronounced oscillation of spin po- 
larization near the potenial boundary x = 0. The peak 
value of Xz is order of 10 -6 meV~ l nm~ 2 , and the peak 
value of P z is about 20%, indicating the SHE here is gi- 
ant. 2). The induced spin polarization is found to be in- 
sensitive to the Fermi energy on the TI surface. This may 
be understood because the spin polarization is approxi- 
mately inversely proportional to Vp, and vf is a constant 
for Dirac particles. This is markedly different from the 
usual 2D system where vf is proportional to \[Ef, hence 
the spin polarization 1/y/Ep. 3). The oscillation period 
at the zero gate region is inversely proportional to the 
Fermi wave vector, or 1/kp = Tivf/Ef, typically tens 
of nm for Ep at tens of meV, and the period at the bot- 
tom gated region is proportional to 7ivf/\Ep — V2I, and 
can be larger. The spin polarization predicted in our the- 
ory may be detected in spin resolved STM on mesoscopic 
samples with over several micron mean free path. Differ- 
ent from the 2D electron gas in semiconductors formed 
in sandwiched interface, which is difficult to use STM, 
the surface of TI can be directly measured by STM. The 
bottom gated device allows to detect the spin polariza- 
tions at both left and right regions. And we notice that 
the back-gated TI thin film device has been realized in 
experiment recently^ 

To further illustrate the giant SHE on the TI surface, 
we compare our result with the theoretically estimated 
SHE in ballistic 2D electron system (2DES) with Rashba 
spin-orbit coupling. Note that the ballistic SHE in 2DES 
is relatively weak and has not been observed in experi- 
ment. In the 2DES case, the Fermi velocity is propor- 
tional to square root of the Fermi energy, so the spin 
polarization is inversely proportional to the square root 
of the Fermi energy^. For a typical semiconductor such 
as InGaAs/InAlAs heterostructure^, we have effective 
mass m* = 0.05m e , with m e the free electron mass, and 
the Rashba spin-orbit coupling an « lOmeV ■ nm . The 
theoretical calculation^ indicates that the peak value of 
P z w 2.4% if E F « 3.3meV, and P z w 0.4% for a more 
realistic Fermi energy Ep « lOOmeV corresponding to 
the 2D electron density of «2D ~ 2.1 x 10 12 cm~ 2 . There- 
fore, the local spin polarization we predicted for the TI 
surface is about 50 times larger than that in a typical 
2D semiconductor with Rashba spin-orbit coupling. We 
also remark that the SHE we predicted in the TI is about 
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FIG. 3: (Color online) Spatial profiles of spin polariza- 
tion density \ z {x;Ef) (top panels), and spin polarization 
P z (x;Ef) (bottom panels) on the surface of TI, described 
by Eq. (1) with up = 5 x 10°m/s. Left panels (a) and (b) are 
for n — n junction: Vi = 0, V2 = 40 meV, and the Fermi en- 
ergy Ef = 50 meV (red dashed line) and Ef = 60 meV (blue 
solid line). Right panels (c) and (d) are for n — p junction: 
Vi = 0, V2 = 60 meV, and -Ef = 35 meV (red dashed lines) 
and Ef = 40 meV (blue solid lines). P z — 1 corresponds to 
completely polarized spins at Ef- 



1000 times larger than that estimated for HgTe quantum 



well where the spin-orbit coupling is induced by in-plane 
potential gradient] 26 ' 34 We conclude that the surface of 
TI should be an excellent candidate to observe the SHE 
in the ballistic limit. 

In summary, we have theoretically examined the meso- 
scopic spin Hall effect on a surface with a step function 
potential in 3-dimensional topological insulator. By ap- 
plying the McMillan Green's function method, which is 
based on the scattering wave function method and was 
previously used in study of superconductor junctions, we 
have derived analytic expressions for the electric current- 
induced spin polarizations on the surface (actually, this 
analytical method is suitable for the ballistic SHE in var- 
ious systems 3 ^,). In the ballistic transport limit, a gi- 
ant spin polarization oscillation across the junction is 
induced by a transverse electric current. The spin po- 
larization is estimated as large as 20%, which is insen- 
sitive to the Fermi level and is not confined by the po- 
tential step due to the Klein paradox. Its magnitude is 
about two orders larger than that in 2-dimensional elec- 
tron gas with Rashba spin-orbit coupling. The spatial 
oscillation period is order of inverse of the Fermi wavevec- 
tor. These features are markedly distinguished from the 
2-dimensional electron gas and should be observable in 
experiments such as spin resolved scanning tunneling mi- 
croscope. 
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